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The physical mechanism responsible for driving accretion flows in astrophysical accretion disks is 
commonly thought to be related to the development of plasma instabilities and turbulence. A key 
question is therefore the determination of consistent equilibrium configurations for accretion-disk 
0^ \ plasmas and investigation of their stability properties. In the case of collisionless plasmas kinetic 

theory provides the appropriate theoretical framework. This paper presents a kinetic description of 
low-frequency and long-wavelength axisymmetric electromagnetic perturbations in non-relativistic, 
strongly-magnetized and gravitationally-bound axisymmetric accretion-disk plasmas in the colli- 
sionless regime. The analysis, carried out within the framework of the Vlasov-Maxwell description, 
relies on stationary kinetic solutions of the Vlasov equation which allow for the simultaneous treat- 
ment of non-uniform fluid fields, stationary accretion flows and temperature anisotropies. It is 
demonstrated that, for such solutions, no axisymmetric unstable perturbations can exist occurring 
on characteristic time and space scales which are long compared with the Larmor gyration time 
J-j ' and radius. Hence, these stationary configurations are actually stable against axisymmetric kinetic 

instabilities of this type. As a fundamental consequence, this rules out the possibility of having the 
axisymmetric magneto-rotational or thermal instabilities to arise in these systems. 
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A fundamental issue in the physics of accretion disks (ADs) concerns the stability of equilibrium or quasi-stationary 
configurations occurring in AD plasmas. The observed transport phenomena giving rise to the accretion flow are 



commonly ascribed to the existence of instabilities and the subsequent development of fluid or MHD turbulence 
OO ■ Grl3- I n principle, these can include both MHD phenomena (such as drift instabilities driven by gradients of the 
t-H ' fluid fields) and kinetic ones (due to velocity-space anisotropies, including, for example, trapped-particle modes, 
cyclotron and Alfven waves, etc.). Possible candidates for the angular momentum transport mechanism are usually 
identified either with the magneto- rotational instability (MPJ) 0-0 or the thermal instability (TMI) f^-fHI > caused by 
unfavorable gradients of rotation/shear and temperature respectively. The validity of the above identifications needs 
t-H . to be checked in this case, because they usually rely on incomplete physical descriptions, which ignore the microscopic 
(kinetic) plasma behavior. In fact, "stand-alone" fluid and MHD approaches which arc not explicitly based on kinetic 
theory and/or do not start from consistent kinetic equilibria, may become inadequate or inapplicable for collisionless 
or weakly-collisional plasmas. Apart from possible gyrokinetic and finite Larmor-radius effects (which are typically 
not included for MRI and TMI) , this concerns consistent treatment of the kinetic constraints which must be imposed 
on the fluid fields (see related discussion in Refs.[l2l. [l3|). This concerns, in particular, the correct determination of 
the constitutive equations for the relevant fluid fields. Because of this, the issue of stability of these systems is in need 
of further study. 

In this regard, some relevant background materials are provided by Refs. [l2H14j |. where a perturbative kinetic theory 
for collisionless plasmas has been developed and the existence of asymptotic kinetic equilibria has been demonstrated 
for axisymmetric magnetized plasmas. In AD plasmas, in particular, they are characterized by the presence of sta- 
tionary azimuthal and poloidal species-dependent flows and can support stationary kinetic dynamo effects, responsible 
for the self-generation of azimuthal and poloidal magnetic fields [15| , together with stationary accretion flows. This 
provides the basis for a systematic stability analysis of such systems. We stress that these features arise as part of 
the kinetic equilibrium solution, and are not dependent on perturbative instabilities. Furthermore, by assumption 
in the theory developed here there is no background (i.e., externally- produced) radiation field. In principle, for a 
collisionless plasma at equilibrium, charged particles can be still subject to EM radiation produced by accelerating 
particles (EM radiation-reaction). However, the effect of these physical mechanisms is negligible for the dynamics of 
non-relativistic plasmas, and therefore they can be safely ignored in the present treatment. 
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In this paper we address the stability of these equilibria with respect to infinitesimal axisymmetric perturbations. We 
restrict attention to the treatment of non-relativistic, strongly-magnetized and gravitationally-bound (see definition 
below) collisionless AD plasmas around compact objects for which the theory developed in Refs.JH, E3] applies. The 
plasmas can be considered quasi-neutral and characterized by a mean-field interaction. Accretion disks fulfilling 
these requirements rely necessarily on kinetic theory in the so-called Vlasov-Maxwell statistical description, which 
represents the fundamental physical approach for these systems. In AD plasmas, electromagnetic (EM) fields can 
be present, which may either be externally produced or self-generated. At equilibrium, they are taken here to be 
axisymmetric and of the general form B< e «) = B^b = B^ q) +B { p q) and E^) = -V$( e<?) (x). Here B^ e?) = Z(x)V<^ 
and Bp'' = V-0(x) x Vip are the toroidal and poloidal components of the magnetic field respectively, with J(x) and 
$(e«) ( x ) being the toroidal current and the electrostatic potential. Furthermore, (R, <p, z) denote a set of cylindrical 
coordinates, with x = (R, z), while (ip, <p, $) is a set of local magnetic coordinates, with ijj being the so-called poloidal 
flux function. The validity of the previous representation for B^ e9 ^ requires the existence of locally nested magnetic 
-0— surfaces, represented by ip = const., while the expressions for ^(x), /(x) and $( e9 ) (x) follow from the stationary 
Maxwell equations. The gravitational field is treated here non-relativistically, by means of the gravitational potential 
<I>G = < J ) g( x )- This means that the electrostatic and gravitational fields are formally replaced by the effective electric 
field Wjf = -V$ e / f , determined in terms of the effective electrostatic potential <5> e J f = <& (e,3 )(x) + ^$ G (x), 
with M s and Z s e denoting the mass and charge, respectively, of the s-species particle (where s can indicate either 
ions or electrons). Based on astronomical observations, the magnetic field magnitudes are expected to range in the 
interval B ~ I0 1 — 10 8 G (l6l - [l8| . This implies that the proton Larmor radius th is in the range 10 -6 — 10 3 cm 
(the lower values corresponding to the lower temperature and the higher magnetic field). Additional important 
physical parameters are related to the species number density and temperature. Astrophysical AD plasmas can 
have a wide range of values for the particle number density n s , depending on the circumstances considered. Here 
we focus on the case of collisionless and non-relativistic AD plasmas assuming values of the number density n s in 
the range n s ~ 10 6 — 10 15 cto~ 3 . For reference, the highest value of this interval corresponds to ion mass density 
Pi ~ 10~ 9 gcTO~ 3 . The choice of this parameter interval lies well in the range of values which can be estimated 
for the so-called radiatively inefficient accretion flows (RIAFs, [HI, For these systems, estimates for species 

temperatures usually lie in the ranges ~ I — lffikeV and T e ~ 1 — lOkeV for ions and electrons respectively. 
Depending on the magnitude of the EM, gravitational and fluid fields, the AD plasmas can sustain a variety of notable 
physical phenomena, the systematic treatment of which requires their classification in terms of suitable dimensionless 
parameters. These are identified with £m,s, £s and o~ s , to be referred to as Larmor-radius, canonical momentum and 

total-energy parameters. Their definitions are respectively: Em,s = rgxf^ ; £ s =~ 



and er, 



Z s e>S> 
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Here, 



r^ s = Vths/fl cs denotes the Larmor radius of the species s, Vths and f2 cs are the species thermal velocity and the Larmor 
frequency respectively, (AL) eq is the characteristic scale-length of the equilibrium fluid fields, v is the single-particle 
velocity and v v = v • i?V</?. Systems satisfying the asymptotic ordering < a s , e s , e, £m,s 1 are referred to as 
strongly-magnetized and gravitationally-bound plasmas [12. 1 1 31] - with the parameters a s ,£ s and em,s to be considered 
as independent while e = max{e s ,s = l,n}. In the following, we shall assume that the poloidal flux is of the form 

if) = i^(x), with ip(x.) ~ 0(e°), while the equilibrium electric field satisfies the constraint t^^T^T^ | ~ O (e). 
This implies that to leading-order $( e9 ' is a function of if) only, while <& e Jf remains generally a function of the 

type <& e /f — $s (Vs 1 ?) ( see Ref.[la|). At equilibrium, by construction, the particle toroidal canonical momentum 
p vs = ^r4'*s = M s Rv v + ^r'lp, the total particle energy E s = Z s e^^ s = H^v 1 ' -\-Z s e<& e J 'f and the magnetic moment 
m' s predicted by gyrokinetic theory are either exact or adiabatic invariants. In particular, the above orderings imply 
the leading-order asymptotic perturbative expansions for the variables and 

^ s =^[l + 0(e s )], (1) 



= [l + 0(a s )], (2) 

while similarly m' s = [1 + (e M , s )]- Fr om here on, we will use the notation that primed quantities are always 

evaluated at the guiding-center. In particular, w' = v — u'h' — ~V e ±f denotes the perpendicular particle velocity in 
the local frame having the effective drift velocity V e ft = -§r x b', while u' = v- b'. In the following we shall also 
assume that the toroidal and poloidal magnetic fields and the species accretion and azimuthal flow velocities scale as 
jg^l ~ O (e) and n^p^r ~ O (e) respectively. 

In validity of the previous assumptions, an explicit asymptotic solution of the Vlasov equation can be obtained for 
the kinetic distribution function fti (KDF). As pointed out in Ref . |13l] . ignoring slow-time dependencies, this is of 
the generic form ft% = /*f (X* s , (4>*s,'&*s))- Here are the invariants X* s = (E s , -0*s,p^ s , m' s ), while the brackets 
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(■0* s , ®*s) denote implicit dependencies for which the perturbative expansions (p} and ([2]) are performed. Therefore, 
fts is by construction an adiabatic invariant, defined on a subset of the phase-space T — x U, with fi C M 3 and 
U = K 3 being, respectively, a bounded subset of the Euclidean configuration space and the velocity space. Hence, /» s 
varies slowly in time on the slow-time-scale (At) eq , i.e. hi/*f ~ j^y-q - In view of the previous orderings holding for 

AD plasmas, this implies also — <C 1, where r co ; jS denotes the Spitzer collision time for the species s. Therefore, 

this requirement is consistent with the assumption of a collisionless plasma. A possible realization of /*f is provided 
by a non- isotropic generalized bi-Maxwellian KDF. As shown in Ref.[13|], f*l determined in this way describes Vlasov- 
Maxwell equilibria characterized by quasi-neutral plasmas which exhibit species-dependent azimuthal and poloidal 
flows as well as temperature and pressure anisotropies. The existence of these equilibria is warranted by the validity of 
suitable kinetic constraints (see the discussion in Ref . [H, EH ) ■ As a consequence, the same equilibria are characterized 
by the presence of fluid fields (number density, flow velocity, pressure tensor, etc.) which are generally non-uniform 
on the ip— surfaces. 

Let us now pose the problem of linear stability for Vlasov-Maxwell equilibria of this type. This can generally be 
set for perturbations of both the EM field and the equilibrium KDF, which exhibit appropriate time and space scales 
{(At) osc , (AL) osc }. Here both are prescribed to have fast time and fast space dependencies with respect to those of 
the equilibrium quantities, in the sense that 

(At)~ ^ (AL)°- 

(At) eq (AL) eq [ h [ j 

with A being a suitable infinitesimal parameter. In the case of strongly-magnetized AD plasmas, to permit a direct 
comparison with the literature, we also assume that these perturbations are non-gyrokinetic. In other words, they are 
characterized by typical wave- frequencies and wave-lengths which are much larger than the Larmor gyration frequency 
fi cs and radius r^ s . This implies that the following inequalities must hold: 

'' L ~ «1, (4) 



(At) osc (AL) C 

with tls = while A must satisfy A 3> <J S , e s , £, £m,s- These will be referred to as low-frequency and long- 

wavelength perturbations with respect to the corresponding Larmor scales. Notice that Eqs.Q and are independent 
and complementary, establishing the upper and lower limits for the range of magnitudes of both (At) and (AL) 
We now determine the generic form of the perturbations as implied by the above assumptions. For this purpose, we 
shall require in the following that the EM field is subject to axisymmetric EM perturbations of the form <5B = V x <5A, 

8E = -V8<f> - i^f, with | S4> (t, | , f ) ,8 A (f , f , {\ j both assumed to be analytic (with respect to ^ and i?) 
and infinitesimal, i.e., such that j ^^) j ; |bT^) | ~ 0(e). This implies that the corresponding perturbations for the 

EM potentials must scale as i^f^i , r^^jj ~ 0(e)0(\), with A^ denoting the equilibrium vector potential. As a 
consequence 



dt 



1 dSA 



-v 



dt c dt 

Similarly, the perturbation of the equilibrium KDF is taken of the general form 



(5) 



Sfs = Sf s ( X* s , (^* s , $* s ), T' T' T ) ' ( ( ») 



AAA 

with y[j ~ 0(e)0(X). It follows that the corresponding KDF (the solution of the Vlasov kinetic equation) must now 
be of the general form 

^ ■& t 



fs = fs \ X mB , [lp* s , j-, -, - ] , (7) 

while, from the Maxwell equations, the perturbations {8(f), 8 A} are necessarily linear functionals of 8f s . However, 
for analytic perturbations of the form ((7]), f s must itself be regarded as an analytic function of rp and $. Therefore, 
invoking Eqs.([T]) and ([2|), the same KDF can always be considered as an asymptotic approximation obtained by 
Taylor expansion of a suitable generalized KDF of the form jj 9e ") = j]f e ' 1 ) (x^, (ip* s , $* s , Y* s ) , jj , with Y* s = 
£^ i n particular, denoting 8f { s 9en} = /i 9en) - f%, it follows that also 8f { s 9en) is such that Sf { s 9en) = 
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6 fa (X*s, (ip*s: &*s, Y rs ) , j) . Then, by Taylor expansion with respect to the variables F* s , the perturbation Sf, 
can be shown to be related to Sf s (defined by Eq.©) by 

Sf ( 9 en) ^ t } f\ eX p-*, 



gen) 



(8) 



where corrections of an d have been neglected. Similarly, invoking again Eqs.([T]) and ©, for the analytic 

perturbations {8<p, SA} we can introduce the corresponding generalized perturbations {54>( gen \ SA^ gen ^. Neglecting 
in the similar way corrections of an< ^ °o"\) ' these are given as follows: 



^* si. I — 



A'A ,exp " 



A' A 



exp 



(9) 
(10) 



Analogous expressions for the corresponding generalized perturbations can be readily obtained. In particular, using 
Eq.©, we get the following representation for 5fs 9en ^: 



Sf 



{gen) _ z7(gen) 



Sf ( s 9en> (X, s ,(^ s ,^ s ,i; s ))exp 11 



(11) 



where, expanding the Fourier coefficient and neglecting again corrections of and q > Sfs 

6f s (^X* s , (ip* s , ^5-, j Y Therefore, in view of Eq.©, for infinitesimal axisymmetric analytical EM perturba- 
tions {S(f>, 5 A}, to leading order in A the Vlasov equation implies the dispersion equation 



0(0 sfigen) 



iuSf s X* s , (V>* s , $* s ), 



A' A 



(12) 



iujq s 



t/j d 



A A 



-ri -r v • 5 A \ 



A' A 



df 



dE s 



Apart from the trivial solution u> = (i.e., a stationary perturbation of the equilibrium), this requires that, for uj ^ 0, 
one must have 



Sf. 



1 



-v-5A 



(eg) 



dE„ 



(13) 



where, by construction, 8f Sl 5(f) and SA are manifestly independent of uj. Hence, Eq. (|13p necessarily holds also when 
\uj\ is arbitrarily small. In this limit |<5</>, 5A, 6f s \ tend necessarily to infinitesimal stationary perturbations of the 

equilibrium solutions. On the other hand, Eqs.©, (|T0|) and (fTTj) show that (Scj), 5A, Sf s \ are always asymptotically 

close to the generalized quantities which are by definition equilibrium perturbations [i.e., 

functions of ^£^±^ 5 ""a*" )]" Since ^he latter again represent an equilibrium and are independent of w, it follows that 

the only admissible solution of the dispersion equation (|13p is clearly independent of uj as well and coincides with the 
null solution, i.e. 



SA 



A' A 
A' A 



■0 $ 

Sf s ( X» s , (^» 8 ,$* a ), -, - 



0. 
0. 
0. 



(14) 
(15) 
(16) 



In summary: no analytic, low- frequency and long-wavelength axisymmetric unstable perturbations can exist in non- 
relativistic strongly-magnetized and gravitationally-bound axisymmetric collisionless AD plasmas. We stress that this 
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result follows from two basic assumptions. The first one is the requirement that the equilibrium magnetic field admits 
locally nested ^—surfaces. The second one is due to the assumed property of AD plasmas to be gravitationally-bound. 
This implies (as pointed out above) that the effective ES potential <&y* is necessarily a function of both ip and 
and therefore the perturbation of the KDF is actually close to a function of the exact and adiabatic invariants X* s . 
A notable aspect of the conclusion is that it applies to collisionless Vlasov-Maxwell equilibria having, in principle, 
arbitrary topology of the magnetic field lines which can belong to either closed or open magnetic ^—surfaces. Also, 
as pointed out in Refs.(l2l. fl3| . for strongly-magnetized plasmas these equilibria can give rise to kinetic dynamo 
effects simultaneously with having accretion flows. These results are important for understanding the phenomenology 
of collisionless AD plasmas of this type. In particular, they completely rule out the possibility that axisymmetric 
perturbations, which are long- wavelength and low- frequency in the sense of the inequalities ([?]), could give rise to 
kinetic instabilities in such systems. This conclusion applies for collisionless AD plasmas (having in particular particle 
densities within the range mentioned earlier) which are strongly-magnetized and simultaneously gravitationally-bound 
. Since fluid descriptions of these plasmas can only be arrived at on the basis of the present Vlasov-Maxwell statistical 
description, also MHD instabilities, such as the axisymmetric MRI @, [13], the axisymmetric TMI (see for example 
0-[ll|), and axisymmetric instabilities driven by temperature anisotropy (e.g., the firehose instability [HI) remain 
definitely forbidden for collisionless plasmas under these conditions. 
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